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Abstract 

We study the low momentum behavior of nucleon gravitational form factors in the framework 
of the heavy baryon chiral perturbation theory. At zero recoil they determine the momentum and 
spin apportion between nucleon constituents. Our result provides an insight into the response 
of the nucleon's pion cloud to an external weak gravitational field and establishes a theoretical 
framework for extrapolation of experimental and lattice data on the nucleon form factors to zero 
momentum transfer. We also discuss form factors corresponding to higher-rank tensor currents 
related to the moments of generalized parton distributions. 
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1. The fundamental problem of hadronic physics is the understanding of the structure of the 
nucleon, — the building block of our planet and its inhabitants, — in terms of quark and gluonic 
degrees of freedom. The nucleon's mass and momentum content is one of the issues being under 
study for decades. Higgs mechanism giving quarks their masses generates only a tiny fraction 
of the total nucleon mass, the bulk of which comes entirely from the complicated dynamics of 
colored glue. This idea was supported by deeply inelastic scattering experiments which revealed 
that quarks and gluons are equally responsible for the energy-momentum structure of the nucleon 
with both of them carrying about 50% of its momentum. A similar study has to be addressed 
for the spin content of hadrons. The nucleon being a composite particle builds its spin from the 
angular moment aQ of its constituents J g>s , i.e., in the case of quarks from their spin AE/2 and 
orbital motion L q , 

1 AS 

2 = — + L i + J 9- (!) 
A next natural question arises about the numerical apportion between these contributions. The 
naive quark model (QM) attributes the whole of the nucleon spin to the sum of the ones of quarks 
ASqm = 1 [0]- However, the experimental data tells the opposite: only a small fraction of the 
total spin is carried by quarks AS cxp = 0.16 ± 0.08 @. Thus in order to test the spin sum rule 
(|I|) one has to have access to the other sources of the nucleon spin L q and J g . The nucleon's 
angular momentum is expressed in QCD by a hadronic matrix element of the coordinate-space 
moment of the quark and gluon momentum densities given by specific components of the energy- 
momentum tensor. At this point recall that the magnetic moment of a particle is the space 
moment of the electric current flow and can be measured as a matrix element of the electric 
current at non-zero momentum transfer, i.e., in elastic electron-nucleon scattering by accessing 
electromagnetic form factors. The same methodology is applicable in the current circumstances 
so that the missing orbital momentum parts can be deduced from off-forward matrix elements 
of the QCD energy-momentum tensor. Thus, in order to understand the nucleon spin content 
one has to study nucleon gravitational form factors. These functions, which arise in the nucleon 

^Tor a gauge particle the decomposition of its angular momentum into its spin and orbital components is 
not possible. The hand-waiving argument goes as follows: In quantum mechanics, the wave function of a spin-s 
particle is a symmetric rank-2s spinor having (2s + 1) components which transform into each other under the 
rotation of the coordinate system. The orbital wave function is related to the coordinate dependence of wave 
functions and is given by the spherical harmonics of order I for the angular momentum I of the system. Therefore, 
in order to distinguish clearly between the spin and orbital momentum, the spin and coordinate properties must be 
independent. As it is obvious, this condition is not fulfilled for gauge particles whose description in terms of field 
operators inevitably involves a gauge condition. For instance, in the Coulomb gauge the gluon wave function is 
given by the three-potential A a (x) equivalent to the second rank spinor which, however, is a subject for the gauge 
condition divA a (ir) = 0. As a result the coordinate dependence of the vector cannot be independently defined for 
each of its components and leads to the inability to separate the spin and orbital degrees of freedom. 
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scattering off a weak gravitational field, are measurable experimentally in exclusive processes 
such as deeply virtual Compton scattering || |], |5| and hard exclusive meson production || [?J . 
Moreover, they can be simulated on a lattice. As we just said, apart from the interest in their 
own right, gravitational form factors at zero recoil determine the quark and gluon contributions 
to the angular momentum. In practical determination of the latter either in laboratory or lattice 
experiments one needs in all cases to perform an extrapolation to the limit of zero momentum 
transfer A 2 = to get the net effect. The experimental kinematics is restricted to low t-channel 
momenta of order 0.1 GeV 2 |J and, thus, favors the applicability of the chiral perturbation theory 
(xPT) || to unravel the small-momentum transfer behavior of the form factors in question. On 
a lattice, the minimal achieved momentum is limited by the size of the lattice. The first Monte 
Carlo simulations ||, 1(| were done with lattices 16 2 x 32 where A 2 ^ ~ 0.4 GeV 2 . Thus in order 



to enter the regime of validity of x?T one needs repeat the simulations with bigger lattices. In the 
present Letter we compute one-loop contributions to off-forward matrix elements of local quark 
and gluon composite operators and thus determine their momentum transfer dependence. 

2. The spin and momentum structure of the nucleon is carried by off- forward matrix elements 
of quark and gluon composite operators. In our consequent presentation we will be concerned 
with the parity-even sector which can shed some light on the spin structure of the nucleon. Under 
chiral SU{2) L ® SU (2) R group the operators transform either as isovector (3, 1) © (1, 3) or isoscalar 
(1, 1) representations and correspond to flavor non-singlet TZ q,a and singlet TZ q, ° quark operators, 
respectively, 

KL-h = w S H $T A j^iV m iV N $ , (2) 

where r A = (1, r a ) with 1 = lpixpi- The operator of symmetrization and trace subtraction is 
defined as follows, e.g., S WAt2 t WM2 = ^(i W(U2 +* jU3W - Ig^^t^). The isoscalar combination Tl q >° 
mixes under the renormalization group evolution with gluonic operators, 



K^-h = ^ S „ M . G^iV^ iV,» . ■ -iV^ G NU . (3) 
The nucleon matrix element of these operators are parametrized via 'form factors' as follows 

(p2\n A \ Pl ) = S u( P2 )r A lfll u(pi) K 2 • ••/y-U.U 2 ) + • • • + A jU2 . . . a^,i(A 2 )} 

where we have introduced the vectors p = P1+P2 arid A = p 2 — p\ and assumed the normalization 
for hadronic states (p\p) = 2poV and bispinors uu = 2M. We have reduced in the above equation 
all tensor structures of the scalar Dirac bilinear containing at least one vector p^ by means of 
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the Gordon identity to the vector and tensor Dirac bilinears of the first two lines. There is an 
important constraints on the matrix elements (f|) imposed by the time-reversal invariance, namely, 
vanishing of Lorentz structures with odd powers of the momentum transfer A ||. We will keep 
this condition in mind without implementing it explicitly in our subsequent considerations. 

The reduced matrix elements A, B and C are related to the moments of parity-even quark or 
gluon generalized parton distributions [JTTJ, 12], |j] measurable in deeply virtual Compton scattering 



via 



/ dxx j ~ 1 H(x,rj, A 2 ) = ^^-^(A 2 ) +7 / >'C j (A 2 ) , 
f dxx^E^x^A 2 ) = Y,V k B J , 1 - k (A 2 )-r ] iC J (A 2 ). (5) 

J —1 7. n 



The angular momentum sum rule (|1|) involves the form factors at zero recoil entering their second 
moment |T^] 

J=^im o {A 2)2 (A 2 )+5 2 , 2 (A 2 )}, (6) 

which are attributed to quarks or gluons depending on the contributing composite operators. 

We consider single-nucleon systems and, in order to have a consistent power counting, use the 
formalism of the heavy baryon x?T which treats the nucleon as a non-relativistic infinitely heavy 
particle ||13|| . For the purpose of a unique translation of the Lorentz-covariant matrix elements in 



Eq. (H) into non-relativistic ones its is convenient to use the Breit reference frame JHL It is 



defined by the condition p = which leads in turn to A = 0, p = 2MJI - A 2 /(AM 2 ) = 2Mb, 
and p 2 = —Pi = A/2. In the Breit frame, the Dirac bilinears arising in Eq. (^) are reduced via 
equations 

u{P2)liJ,u{pi) = v II u v (p 2 )u v (p 1 ) + ^u v (p 2 )[S^ S ■ A]u v {p\) , 

u(P2) Za ^ V u(pi) = V^^^Uv^Uvip!) + ^u v {p 2 )[S^ S ■ A]u v ( Pl ) , 

to the ones constructed from the projected large components of the heavy baryon bispinor u v (pi) = 
(1 + yf)u(pi) I \j2-\- 2{y ■ pi) /M normalized as u v u v = 2M. The velocity v arises in the heavy-mass 
decomposition of the incoming and outgoing nucleon momenta pi t2 = Mv +fc=pA/2 and the residual 
momentum is defined as k = (M(S — 1),0) when v = (1,0). The spin vector = §0^1,751^ 
reduces to the spin three- vector S — (0, 27/2) in the nucleon rest frame. Exploiting these results, 
the decomposition @ transforms in the Breit frame to the equation 



(Pal^l. Jpi> = S«,b 2 )rVbi)^ 1 {(2M5)^ 1 i? JJ (A> )12 . . . v N + ... + %i(A 2 )A M2 . . . A Mj } 
+ S ±-u v {p 2 )r A [S^S ■ A]u v ( Pl ) {(2M5y- l M hJ (A 2 )v^ ...v H + ... + M ]A (A 2 )A, 2 ...A N } 
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+ — — A * • • • A « Cj ' ' (7) 

where we have introduced a short-hand notation for the generalized electric- and magnetic-like 
form factors 



E 3,k = A j,k + -7TH B j,k , M jyk = A jtk + B j>k 



A 2 

AM 2 ' 

Now we are in a position to address the calculation of the small-A dependence of these form factors. 
As it is obvious from the tensor decomposition, one cannot compute form factors accompanying 
Lorentz tensors involving several A M , i.e., k < j — 2, while restricting oneself to next-to-leading 
order in the chiral expansion. Also note that E\\ and correspond to the standard Sachs 
nucleon electromagnetic form factors extensively discussed in the literature |14], pl| and, therefore, 
their evaluation will not be repeated presently. 

3. The calculational procedure is rather straightforward. First, one constructs composite 
operators from the nucleon and pion fields which match the quantum numbers of the ones on 
the quark-gluon level and adds them to the effective Lagrangian [0, [I7| . Then one uses 



the heavy-baryon chiral perturbation theory fL3| , see also the review |18j, for the computation of 
tree and one-loop contributions. For the case at hand, one needs the leading order pion-nucleon 
Lagrangian 

f 2 

C = N v {iv ■ V + 2g A S ■ A} N v + ^tr {d„S d^) + A tr [M q [E + E ] ) } . (9) 
The axial-vector and vector (in the covariant derivative = + V^) pion potentials are 

= i (?d^ + £d^) , A. = \ (? d ^ - ^) ■ ( 10 ) 

Here E is constructed from the multiplet of the pion fields 

E = e = exp (jn ■ r) , (11) 



and f n = 93 MeV the pion decay constant. Expanding to the first non-trivial order in the 
pion fields we get for the axial and vector potentials: = —-^-8^ ■ r + ... and V M = 
£abc TC 'x a du7r b + . . ., respectively. In the mass term, the coefficient A is related to the quark 



condensate via A = —^(ipip) an d the quark mass matrix is M q = diag(m n , m^). 

Ultraviolet divergences generated by one-loop diagrams have to be absorbed into coefficients 
of next-to-leading order counterterms in the chiral expansion. 

Now we are in a position to discuss the construction of a basis of hadronic twist-two operators 
corresponding to the quark and gluon ones in Eqs. (0) and (|3]). There are two types of local 
hadronic operators: built either from the pions alone or bilinear in the nucleon fields. The nucleon 



4 



operators have the form which mimics the tensor decomposition of the off-forward matrix elements. 
The leading operators contributing to isoscalar and isovector combinations of Ej j and Mjj are 

= <(2M)- 1 V . . . v N N v rf + N v 

+ bf(2My-\-id u ) « w . . . v^NvTfJt^N,, + ..., (12) 

where 

T& = ~ (£t a ? + ^) . (13) 

The matching coefficients a N and b N from the partonic to the hadronic level are unknown and have 
to be determined from experimental data. Subleading operators are deduced from this expression 
by a mere replacement of vs in the Lorentz structure by the derivatives. However, due to the 
time-reversal symmetry restrictions one has to replace an even number of velocities. Here and 
everywhere later, the nucleon operators are normalized as follows N v (0)\p) = u v (p)\0). 

The pion sector has to be considered separately for the isoscalar and isovector components. As 
we will see momentarily, one can essentially reduce the number of operators at each level in the 
derivative expansion provided one uses all underlying symmetries. A number of relations arises 
on the basis of the following 'magic' property of the Pauli matrices r 2 r r 2 = — f T which results 
into a useful conjugation property for the nonlinear pion field £, 

t 2 EW 2 = E t . (14) 

4. Let us start from the isoscalar sector which transforms as (1,1) under chiral SU(2)l <8> 
SU(2)r- The leading nucleon operators are displayed in Eq. ( |l"2D for A = so that t° + = 1. Next, 
one can immediately conclude that the lowest local operator in Eq. (|2]) does not correspond to 
any pion operator as a consequence of the unitarity of £ and Eq. flT4|), 

tr{S d^) = 0. (15) 

Its generalization for an arbitrary odd j = 2k + 1 reads 

tr K • • • d ^ - d H S *} = \ d^tr {d, 2 . . . d, k+1 S d, k+2 . . . d^tf } ■ 

(16) 

Therefore, all pion operators with an odd number of derivatives j = 2k + 1 can be reduced to the 
ones containing an odd number of total derivatives, so that we have 

'k-l 



S . \^hi + A-id) 2l+X tx{{-idf- l E (id) k - l tf } 



. 1=0 



fil—fij 



Here, on the right-hand side of the equation, we imply that the Lorentz indices are attached to 
the partial derivatives d in a particular order, e.g., {5 3 } MlAt2M3 = d^d^d^. Note, however, that 
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matrix elements of these operators violate the time reversal properties of the off-forward matrix 
elements in question @ and thus cannot contribute. 

Using the symmetry properties established above we can easily write composite pion operators 
for even j = 2k, 

Ol ni2 ... N = fl$^ {l>- 2 , {-id) 21 tr{{-id) k - l £ {id) k - l Z^ } j . (17) 

Since quantum numbers of hadronic operators match both quark and gluon ones in Eqs. (fj) and 
(|3]), the difference between the latter on the hadronic level arises only in the value of the matching 
coefficients which appear in two species, ajf and aj'f . Below, the index which attributes hadronic 
operators either to the quark or gluon sector will not be displayed explicitly for brevity. Note 
that aj = aj survive in the forward matrix elements and, thus, are measurable in conventional 
deeply inelastic scattering experiments. For instance, a\ is related to the momentum fraction of 
the pion carried by quarks and gluons forming it. For an on-shell pion at low normalization point 
ji 2 ~ 0.3 GeV 2 one has |19| a 2 ' 9 = (x q ) w ~ 0.7 and respectively a^' 9 = (x g ) w ~ 0.3 due to the 



momentum conservation. There is a number of pion operators with more than two E fields since 
their insertion does not alter the twist of the composite operator, however, they do not contribute 
to the one-loop matrix elements we are interested in and will be totally omitted. 

To summarize, the isoscalar nucleon gravitational form factors, — the second moment of the 
generalized parton distributions H and E, — receive contributions from the following pion and 
nucleon operators: 

<3U = flal S tr {d^E } , (18) 

0^ 2 =a»(2M) S v„ lV , 2 N v N v + b»(2M)(-id v ) S v^J^^-N, 

H'lH'2 ^1/^2 jyl 

+ Cil^Sj-^)(-^)NvN v , (19) 

and we have also included an operator with two total derivatives which generates the structure 
C 2 (A 2 ). 

A calculation of the diagrams (a) and (c), (b) and (d) do not contribute for isoscalar operators, 
gives 

£ 2 , 2 (A 2 ) = < + { ~ ^) Jo dX + t^} ' 

M 2 , 2 (A 2 ) = b» + 3^j 2 jf 1 dx (a> 2 (x) - # m 2 5(x)) jln ^ - l| , 

C 2 (A 2 ) =c» + Zal-f-M (2m 2 - A 2 ) C dx , (20) 

167r / 2 v J J o Jm 2 (x) 
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(b) 



(d) 



Figure 1: One loop diagrams contributing to the matrix elements of the twist-two operators. The 
self-energy insertions into the external line, which are not displayed explicitly, have to be added. 

where m 2 (x) = m 2 — x(l — x)A 2 . Note the absence of one-loop contributions to the E and 
C form factors from the nucleon operators. It is a mere consequence of a cancellation of the 
one-loop contribution (c) by self-energy insertions into the external lines. We have minimally 
absorbed ultraviolet divergences generated in M22 into higher order counterterms that results 
into a redefinition of their coefficients, with which they enter in the effective Lagrangian —>■ 
rj + coeff L e , with L £ = ~ — 7^ + ln47r. So that the right-hand side of M 2j2 has extra analytic 
terms 

m 2 A 2 

r r i + 7Z — T^r 2 ■ 



(AnU) 2 (4vr/ OT 

Note, that one can also completely absorb analytical contributions present in ( PU[ ) into coun- 
terterms. This simply sets the renormalization scheme prescription which has to be used for all 
observables once are fitted to experimental data in one. 

Another comment concerns E 2j2 - Pion operators generate 1/M-suppressed contributions to 
the structure of interest. Therefore, by power counting we have to add 0(1/ M) bilinear nucleon 
operators constructed from one large and one small components of the nucleon field as well as 
analogous terms stemming from the chiral Lagrangian. Lorentz invariance fixes unambiguously 
their coefficients and no new low energy constants arise. However, the effect of these contributions 
in form factors is A 2 -independent and analytic in the pion mass. Thus they will not computed by 
us presently, although can be anticipated to generate an addendum — 5a 2 gA m3 / (327r/ 2 M) on the 
right-hand side of £"2,2 i n Eq. (|20|) required by momentum sum rules discussed below. Also, there 
arises a Foldy-like term (b^ — a^)A 2 /(4M 2 ) in Ejj, similarly to the form factor case [15], from 
contributions of the small nucleon field components in the heavy-mass expansion of relativistic 
nucleon operators. 

Sum rules for the total momentum and spin of the nucleon (pion) impose constraints on the 
coefficients a 2 and b 2 {a 2 )-, 

<%« + = 1 , + b^ = 1 , d*« + = 1 . (21) 
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The latter two equations imply that the total gravitomagnetic moment of the nucleon vanishes, 
i.e., S« 2 (0) + Sf i2 (0)=0 

The leading (j, j)-structures of the higher j-moments (j > 2) do not receive non-analytic 
contributions in the momentum transfer at next-to-leading order in the chiral expansion due to 
the absence of relevant pion operators. Thus one gets 

^(A 2 ) = < + -.., M,,(A 2 ) = bf (l - 3^ In ^ + . . .) , (22) 

where ellipsis stand for, at least, m 2 /(47r/ 7r ) 2 -suppressed analytic contributions from one-loop 
diagrams and counterterms. 

5. In the isovector sector, the use of Eq. ( |HD results into 

JSte^d^E d^}=0. (23) 

Its generalizations for even j = 2k are straightforward and read 

^i X {^...d llh Ed ltk+l ...d lli ^} = Q. (24) 

Using this property it is easy to convince oneself that pion operators for even j = 2k are reduced 
to the ones involving an odd number of total derivatives, translated in the momentum space to 
odd powers of the momentum transfer, 

m S M . {E4 ;0) H<9f +1 trr" {{-idf- l Z {idf- l+1 Z^ + {-id) k ~ l tf {id) k - l+1 z}\ 



. 1=0 



which, therefore, cannot contribute to the matrix elements in question due to the time-reversal 
condition alluded to above. 

For odd j — 2k + 1, an elementary consideration along the same line as above leads to the 
following structure of the non-singlet pion operators 

<«.-* = f ,, s w ' 25 > 



x 



k 

E4 ; S (-id) 2l trr a {{-id) k ~ l E (id) k ' l+1 E^ + {-idf- l E^ (id) k - l+l E) 

.1=0 



111... II j 



The leading nucleon operators are given in Eq. flEp with A = a. 

Thus, for the isovector gravitational form factors there no contributions from coupling to pions 
while the nucleon operators read 

OX = a»(2M) v^Nr^N + b»(2M)(-id u ) v^J^^-N 
1 

2M 



+ C 2WT7 S i-id^i-id^Nr^N , (26) 
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with r ? a + = \ (fr a ft + ^r a ^). 

Due to the absence of the pion cloud contribution at this order of x?T, i.e., diagrams (a) and 
(b), no non-analytic dependence on the momentum transfer arises. However, the diagrams (c) 
and (d) develop chiral logarithms in the pion mass of the form 

f Tfl 2 I TYl \ 1 

_^ + Dln I| -^-lj}, 

Tfl 2 I 7T1 2 \ I 

f (3gl + l)ln^-gl-l)\ . (27) 



£ 2 , 2 (A 2 ) 


= a N 2 


M 2 , 2 (A 2 ) 




C 2 (A 2 ) 


= c N 2 



K (47T/ W ) a V A X 

We imply that one adds counterterms to the right-hand side of these equations linear in m 2 whose 
(unknown) coefficients absorb minimally the ultraviolet divergences stemming from loops, coeff L e , 
and the Foldy-like term as discussed after Eq. (^). The leading structures of the higher moments 
Mjj and Ejj, apart from the change of an overall normalization b 2 — > b^ and a 2 — > of , have the 
same dependence on the chiral logarithms as in Eq. (|27p. 

6. In this Letter we used the framework of the heavy baryon chiral perturbation theory in one- 
loop approximation in order to predict the small momentum transfer dependence of the moments 
of generalized parton distributions. At next-to-leading order in the chiral expansion only the first 
two moments develop a non-analytic A 2 behavior: isovector electromagnetic form factors and 
isosinglet gravitational form factors. Higher moments can have at most a linear A 2 -dependence 
at this order of x?T. 

Although, the momentum and spin sum rules suggest that the total quark and gluon anomalous 



gravitomagnetic moment is zero p(| ^TJ, the gravitomagnetic radius of the nucleon does not vanish 
and reads 

(r 2 ) = 6 



rfM 2 , 2 (A 2 ^ ~ 2 f - 2 



_- 3a; *ri +1 } +6 (w (28) 



dA 2 

A pragmatic application of our present results consists of their use for extrapolation of experi- 
mental measurements of the angular momentum sum rule (|B]) to zero momentum transfer || and 
as a fit formula for lattice simulations of the parton angular momentum by means of the nucleon 
gravitational form factors 0. Our considerations can be extended to other single- nucleon observ- 
ables such as parity-odd generalized parton distributions as well as off-forward quark transversity 
and tensor gluon form factors. 

We would like to thank J.-W. Chen and T.D. Cohen for useful discussions. This work was 
supported by the US Department of Energy under contract DE-FG02-93ER40762. 
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